Critical behavior in vacuum gravitational collapse in 4 + 1 dimensions 
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We show that the 4 + 1 dimensional vacuum Einstein equations admit gravitational waves with 
radial symmetry. The dynamical degrees ol freedom correspond to deformations of the three-sphere 
orthogonal to the (t, r) plane. Gravitational collapse of such waves is studied numerically and shown 
to exhibit discretely self-similar Type II critical behavior at the threshold of black hole formation. 



INTRODUCTION AND SETUP 



The gravitational collapse to a black hole is a sub- 
ject of intensive studies in general relativity. One of the 
main goals of these studies is to prove the cosmic cen- 
sorship conjecture which says that a physically realistic 
generic gravitational collapse cannot result in a naked 
singularity. It would be most interesting to assert that 
this conjecture is true in vacuum, however with current 
analytical techniques the problem seems tractable only in 
spherical symmetry and in this case there is no vacuum 
collapse because of Birkhoff 's theorem. Thus, in order 
to generate spherically symmetric dynamics one has to 
couple matter fields. A simple choice, which led to im- 
portant insights, is a real massless scalar field. For this 
matter model Christodoulou showed that for small ini- 
tial data the fields disperse to infinity [l| , while for large 
data black holes are formed Q. The transition between 
these two scenarios was explored numerically by Chop- 
tuik leading to the discovery of critical phenomena at 
the threshold of black hole formation. Similar phenom- 
ena were later observed in many other matter models in 
spherical symmetry (see |J| for a comprehensive review) 
but, because of numerical difficulties, only once in vac- 
uum for axially symmetric gravitational waves Q • 

The aim of this letter is to show that - at the price 
of going to 4 + 1 dimensions - one can evade Birkhoff 's 
theorem and have gravitational collapse of pure gravita- 
tional waves in radial symmetry. The idea is very sim- 
ple and is based on the fact that the geometry of the 
three-sphere S 3 has the property that one can break the 
isotropy but still have a homogeneous space. This hap- 
pens as follows. The group of rotations acting on S 3 in 
Euclidean space has a subgroup G3 acting simply transi- 
tively on the three-sphere. This subgroup is isomorphic 
to the universal covering group of the connected com- 
ponent of the rotation group in three dimensions. The 
action of G3 is generated by the simultaneous rotations 
in the (x — y,z — w)-planes, the (x — z,y — u>)-planes and 



(x — w, z — y)-planes (this action on S 3 defines also the 
Bianchi IX homogeneous cosmological model in the 3 + 1 
dimensional general relativity). Thus, in 4+1 dimensions 
it is consistent to consider spacetimes with the metric of 
the form 

3 

ds 2 = -U(t,r)dt 2 +V(t,r)dr 2 +Y,L 2 k (t,r)<jl (1) 

where a k are three one-forms invariant under G3 satis- 
fying the commutation relations d<7i — &j A a k - In 
terms of Euler angles (0 < 6 < w, < (f>, ij) < 2tt) 

ai +icr 2 = e^(cos0 d(j>+id9), 03 = dip - sin d(f>. (2) 

The metric functions Lk(t, r) are the three principal cur- 
vature radii of the squashed three-sphere. The case 
when all three L k are equal corresponds to the standard 
spherically symmetric ansatz for which the Birkhoff the- 
orem applies and the only solutions are Minkowski and 
Schwarzschild. However, if L k are different we will ob- 
tain nontrivial vacuum solutions with gravitational ra- 
diation. In this letter we restrict ourselves to a special 
case of the ansatz in which L\ = L2. Using the 
coordinate freedom in the two-space orthogonal to the 
group orbit of G3, we choose the volume radial coordi- 
nate r = (f;oZ(iS ,3 )/27r 2 ) 1 / 3 , and write the metric as 

)r 2 [e 2B (a 2 



ds 2 = -Ae- 2S dt z +A- l dr z 



'2) ' u 3i j 

(3) 

where A, S, and B are functions of t and r. Note that 
for this ansatz the three-sphere has a residual isotropy of 
the twisted product S 2 x S 1 (as in the Taub universe). 

Substituting the ansatz into the vacuum Einstein 
equations we get the equations of motion for the func- 
tions A(t,r),S(t,r) and B(t,r). In the following we use 
overdots and primes to denote dt and d r , respectively. 
The hamiltonian and momentum constraints are 

' 1 — (8e- 2B ~2e- SB )~2r{e 2S A- 1 B 2 + AB' 2 ^) 

-ArABB'. (5) 
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The evolution equation for B has a form of the quasilinear 
wave equation 

(W~V 3 i?) '-( e ~ s Ar 3 B')'+-e- s r (e~ 2B - e^ B ) = 0. 

(6) 

In addition, we have a slicing condition for 5 

5' = -2r (e 2S A- 2 B 2 + B' 2 ) . (7) 

It is clear from the above equations that the only dy- 
namical degree of freedom is the field B which plays a 
role similar to the spherically symmetric scalar field in 
3+1 dimensions. If B = 0, it is easy to verify (Birkhoff's 
theorem) that the only solution is Schwarzschild 5q = 0, 
Aq = 1 — r 2 /r 2 (or Minkowski if r/j = 0). As we shall 
see below these two static solutions play the role of at- 
tractors. Note that equations (4-7) are scale invariant 
which excludes existence of regular asymptotically flat 
static solutions. 

It is convenient to introduce the mass function m(t, r) 
defined by A = 1 — m(t,r)/r 2 . Then, the hamiltonian 
constraint (4) takes a simple form 

ml = 2r 3 (e 2S A- l B 2 + AB' 2 ^j+-r (3 + - 4 e - 2B ) . 

(8) 

Note that the right hand side of equation © is mani- 
festly positive so m(t, r) is monotone increasing with r. 
For asymptotically flat spacetimes TOqo = lim r ^oo m(t, r) 
exists and is constant in time. The total mass is given 
by M = (STr/SGOmoo. 

We consider the initial value problem for the above 
equations. To ensure regularity at the center we impose 
the boundary conditions 

B(t,r) ~ b(t)r 2 , 1 - A(t, r) ~ a(£> 4 . (9) 

We normalize time by the condition S(t,Q) — 0, which 
means that t is the proper time at the center. 



gaussian 

B(0,r)=p (j-^j e- {r ' ro ^' s \ P(0,r)=rB'(0,r)/r , 

(10) 

where the parameter p was varied and the parameters tq 
and s were fixed. To check universality of the critical be- 
havior, we have verified that the results are independent 
of the specific choice of initial data. 

For all families of initial data the same picture, similar 
to the massless scalar field collapse, emerges. For small 
values of the control parameter p the fields disperse leav- 
ing behind the flat spacetime. For large initial data a 
Schwarzschild black hole is formed, as shown in Figs. 1 
and 2. 









f 




: 








t=0.0 




- t=4.2210314505 




t=4.2703837020 




t=4.2703988129 


- If />* 


- t=4.27040 16743 







0123456789 



FIG. 1: Formation of a black hole for highly supercritical 
initial data. We plot the mass function m(t, r) at the initial 
time and at four late times. The total mass (measured in 
units 3n/8G) is 0.67 (see the plateau of the initial profile). 
During the evolution the function m(t,r) develops a second 
inner plateau which indicates formation of the Schwarzschild 
black hole with mass Mbh = 0.33. 



NUMERICAL RESULTS 

We have solved the above equations using the free 
evolution scheme in which A(t, r) is updated using the 
momentum constraint (5). The hamiltonian constraint 
(4) was only monitored to check the accuracy of the 
code. The wave equation (6) was rewritten as the pair of 
two first order equations for B and an auxiliary variable 
P = e s B/A. Integration in time was done by a mod- 
ified predictor-corrector McCormack method on a uni- 
form spatial grid. The ordinary differential equation @ 
was solved with the fourth order Runge-Kutta method. 
The whole procedure has accuracy of second order in time 
and fourth order in space. 

The numerical results presented below were produced 
for initial data of the form of an 'ingoing' generalized 
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FIG. 2: The plot of B(t,r) for the same data as in Fig.l. 
Outside the horizon developing at r h = 0.57, B tends to zero. 
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The process of settling down to the Schwarzschild black 
hole can be described in more detail using the linear per- 
turbation theory. Linearizing equations (4-7) around the 
Schwarzschild solution we obtain the linear wave equa- 
tion for the perturbation 6B(t,r) 



5B-^rA Q (r 3 A Q SB' 



8A 



SB = 0, A = 1 - 



1 



(11) 



where we have used the scaling freedom to set the radius 
of the horizon = 1. Introducing the tortoise coordinate 
x = r + ~ In £=1 and substituting 5B(t, r) = e~ tkt u(x) 
into we get the Schrodinger equation on the real line 

— oo < x < oo 



where 



-prr + VMx))u = k 2 u, 
ax z 



V(r) - 1 ( 1 M ( 35 + 9 



(12) 



(13) 



Quasinormal modes (i.e. solutions of equation l|12f> sat- 
isfying the outgoing wave conditions u ~ e ± lkx for x — > 
±oo) for potentials of this type have been computed in 
7] via the method of continued fractions. The potential 
()13|l corresponds to the gravitational tensor perturbation 
with I = 2 and in this case the least damped mode (see 
Table III in Q) has eigenvalue k = 1.51 - 0.36 i ( in units 
TZ )■ This mode is expected to dominate the interme- 
diate asymptotics of local convergence to Schwarzschild. 
In Fig. 3 we show evidence confirming this expectation. 



We turn now to the description of critical behavior at 
the threshold of black hole formation. In by now routine 
procedure for bistable systems, using bisection we have 
determined a critical parameter value p* which separates 
black hole formation from dispersion. The behavior of 
near-critical solutions clearly indicates the existence of 
a discretely self-similar critical solution with the echoing 
period A ~ 0.47. The profile of a near-critical solution 
and evidence for discrete self-similarity are shown in Figs. 
4 and 5. 
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FIG. 4: The late time profile of B for a near-critical solution. 
We easily get many echoes even for a moderate fine-tuning 
because the product of the eigenvalue A of the growing mode 
about the critical solution and the echoing period A is quite 
small AA ~ 2.86 (for comparison, this product is about three 
times greater for the massless scalar field critical collapse 3] ) . 
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FIG. 3: Quasinormal ringing of the Schwarzschild black hole. 
We plot the time series In \B(t, ro)\ at ro = 5 for the same 
data as in Fig. 2 but using the proper time at infinity. The 
horizon forms at — 0.575, hence according to the least 
damped quasinormal mode has eigenvalue k — 2.62 — 0.62 i. 
The linear perturbation regime sets in at t ~ 11. The fit of 
an exponentially damped sinusoid to the data on the time 
interval 13 < t < 18 gives k — 2.56 — 0.61 i, which is in good 
agreement with the theoretical prediction. 



0.2 



0.15 



0.1 



m 0.05 



-0.05 



-0.1 
0.001 



A ~ 0.47 




B(ln(r), T q ) 
B ln(r +A, T 
Bln(r ' 
B(ln(r) + 3 A, T 3 ) 



0.01 



0.1 

r 



10 



FIG. 5: Evidence for discrete self-similarity of the critical 
solution. For a near-critical evolution we plot the profile of 
B at some arbitrary central proper time To and superimpose 
the next three echoes which subsequently develop. The times 
T\ , T2 , T3 and the echoing period A were chosen to minimize 
B(ln(r) + nA, T n ) - B(ln(r), T ). 
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As expected, the mass of the black hole Mbh(p) 
changes continuously with p and tends to zero for p — > p* 
according to the power law 



M BH ~ (p* - p) • 



(14) 



where the scaling exponent 7 ~ 0.3289 is universal (i.e. 
independent of initial data). This is shown in Figs. 6 
and 7. Note that, in four space dimensions, mass has 
the dimension of length 2 , hence 7 = 2/A, where A is an 
eigenvalue of the growing mode of the critical solution. 
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FIG. 6: Black hole mass scaling. For supercritical solutions 
the logarithm of black hole mass Mbb (measured in units 
37r/8G) is plotted versus the logarithmic distance to critical- 
ity. The slope of the linear fit yields 7 w 0.3289. The wiggles, 
which are imprints of discrete self-similarity, are shown in 
more detail in Fig. 7. 
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FIG. 7: Fine structure of black hole mass scaling. The linear 
fit from Fig. 6 is subtracted from the data. The period of 
the wiggles agrees to two decimal places with the theoretical 
prediction A/7 = 1.43. 



Summarizing, we have studied gravitational collapse of 
pure gravitational waves in 4 + 1 dimensions and have 
found strong evidence for the Type II discretely self- 
similar critical behavior at the threshold of black hole 
formation. As far as we know, besides the notable paper 
H, this is the only example of critical behavior without 
matter. 

To conclude, let us mention some natural extensions of 
the study presented here. One interesting possibility is 
to investigate the general ansatz Q with two dynamical 
degrees of freedom - the studies in this direction are in 
progress and will be reported elsewhere. It is also nat- 
ural to look for similar models in higher dimensions. If 
one insists that the subgroup of the orthogonal group 
acts simply transitively, the only possibility is the one 
discussed in our paper. However, if one allows multiply 
transitive subgroups of the orthogonal groups, the results 
given in Besse show that there are models similar to 
the one considered here on any odd-dimensional sphere. 
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